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Abstract. In the present we introduce a concept of doubly stochastic 
quadratic operator. We prove necessary and sufficient conditions for 
doubly stochasticity of operator. Besides, we prove that the set of all 
doubly stochastic operators forms convex polytope. Finally, we study 
analogue of Birkhoff 's theorem for the class of doubly stochastic opera- 
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1. Introduction 
Throughout the paper we will consider the simplex 

m 

S" 1 - 1 = {x = Oi,x 2 ,- • • , x m ) € R m : Xi > 0, Vi = 1~m ^x, = 1}. 

i=i 

For any x = (x±,X2 ■ ■ ■ x m ) € S m ~ 1 due to [9], we define x\ = (x[i], x^}, • • • £[m])> 
here xm > xm > ■ ■ ■ > £[m] _ nonincreasing rearrangement of x The point x± 
is called rearrangement of x by nonincreasing. Recall that for two elements 
x, y taken from the simplex S m ~ l we say that an element x majorized by y 
(or y majorates x), and write x -< y(or y >- x) if the following hold 

k k 

^2 x i { \ - ^2 ywf° r an y k = 1 , m ~ 1 - 

i=i i=i 

Note that such a term and notation was introduced by Hardy Littlewood 
and Polya in [7]. It is easy to see that for any x £ 5 m_1 we have 

(-,-,••■ (i,(V--o). 

mm m 

A matrix P = (j>ij)i j=Tm ^ s cane d doubly stochastic (sometimes bistochas- 
tic), if 

Pij > 0, Vi,j = l,m 

m m 

^Pij = 1, Vj = T7m ^Pij = 1, Vi = l,m. 
«=i i=i 
It is known [9] that doubly stochasticity of a matrix P is equivalent to 
Px -< x for all x € S*™" 1 . It is clear that the set of all doubly stochastic 
matrices is convex and compact. Therefore, it was a problem due to Birkhoff 

l 
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concerning description of the set of all extremal points of such a set. A 
solution of that problem was given in [2], and it states that the extremal 
points consist of only permutations matrices. 

The paper is devoted to the same problem mentioned above, but in a class 
of nonlinear operators. To state the problem let us recall some notions. 

An operator V : S" 1 ' 1 — > S m ~ l given by 

m 

{Vx) k = ^2 Pij,kXiXj. (1) 

is called quadratic stochastic operator (q.s.o.), here the coefficients Vij.k sat- 
isfy the following conditions 

m 

Pij,k = Pjt,k > o, y]pzj,fc = l. 

k=l 

One can see that q.s.o. is well defined, i.e. it maps the simplex into itself. 
Note that such operators arise in many models of physics , biology and so 
on. A lot of papers were devoted to investigations of such operators (see, 
for example [1, 3, 8, 10, 11, 12]). We mention that in those papers authors 
studied a central problem in the theory of q.s.o., namely limit behavior of 
the trajectory of q.s.o. 

We say that a q.s.o. (1) is called doubly stochastic (here we saved the 
same terminology as above) if 

Vx^x (2) 

for all x e S™- 1 . 1 

The main object of the paper is the set of all doubly stochastic qua- 
dratic operators (d.s.q.o.). The goal is to study necessary and sufficient 
conditions for doubly stochasticity of quadratic stochastic operators and ex- 
tremal points of such the set of all doubly stochastic quadratic operators. 
To do it, first in sec. 2 we study some properties of doubly stochastic q.s.o. 
Then we give necessary conditions for doubly stochasticity. In section 3 we, 
somehow, describe the set of doubly doubly stochastic operators. Finally, 
in section 4 first of all we study a sufficient condition for q.s.o. to be doubly 
stochastic and then study extreme points of the set of such operators. Note 
that a part of the results was announced in [4], [14]. 

2. Necessary conditions for doubly stochasticity of operators 

This section is devoted to some properties of the doubly stochastic q.s.o. 

Theorem 2.1. Let V : S"* -1 -► S" 71 " 1 - be a d.s.q.o.. Then the coefficients 
Pij t k satisfy the following conditions 



Note that in mathematical economics such an operator is usually called the operator 
of prosperity. 
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m 

a ) ^2 Pi i' k = m > Vfc = m - ( 3 ) 

m 

b ) J2 Pi ^ k -2 } Vi ' fc = V™- ( 4 ) 
here I = {1, 2, • • ■ , m}, |a|— cardinality of a. 

Proof, a) Let C = ^, • ■ ■ , ^). From the definition we have F(C) -< C, 
on the other hand C -< x, Vx € S" 11 ^ 1 . Therefore V(C)j, = C| and since 
Cj, = C one gets V(C) = C which implies (3). 

c) Let Vx -< x. Then according to [Ma] there exists a doubly stochastic 
matrix P(x) = (pij(x)) (which depends on x) such that Vx = P(x)x. Put 
x°(a) = {x\, x®, • • • , x^J, nere 

x i = h>* G a 
x® = 0,i a, 

where a is an arbitrary subset of /. 

It is evident that \a\ < m, therefore x° G S""" 1 and we have 

m ^ \a\ m ^ ^ 

i,j£a i=l i=l 

Since a is an arbitrary set ,then we infer that 

^ Pij,fc < |a|, Va C /, VA; = l,m. 

Using c) we are going to prove b). 

TO 

b) We will prove Yl Pioj k > \ for some fixed io- 
From the equality (see (3)) 

m m m 

m =J2Pij,k= ^ Pij,k + Pi io,k + 2 ^ Pw'.* 

and 

TO 

^ Pij,fc < m - 1 

we obtain 

TO 

or 
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m 

Therefore, b) is satisfied. □ 

3. Description of the class of doubly stochastic operators 

In this section we are going to describe a class of doubly stochastic oper- 
ators. 

Recall that a matrix T = (%), i, j = 1, m is said to be stochastic if Uj > 

m 

and Yl f ij = 1- 

j=i 

Let ^4 = (cijj), i,j = l,m be a symmetric matrix, with > 0. 
Consider the following equation with respect to T: 

A=i(X + n (6) 

here T" is the transposed matrix. 

Below we are going to study conditions for solvability of equation (6) in 
the class of all stochastic matrices. 

Let G m be the group of permutations of m elements. For g € G rn by A g 
we denote a matrix A g = (a g (i) g (j)), hj = l,m, which is called row and 
column permutation of A. 

The following assertions are evident: 



(Ag) g -1 = A 

for any g £ G m 

ii) If A symmetric, then A g is also symmetric for any g € G m 
in) If A stochastic, then A g is also stochastic for any g G G m 
iv) 

(A + B) g = A g + B g ; (XA)g = XA g . 

From i) — iv) we conclude that if T is a solution of (6), then T g is a 
solution of the equation 

A g = \{T + T') 
At first we will study the following set 

Ui = {A = (a^) : a^ = aji > 0, ^ <\a\, ^ = m}. 

i,jea i,jei 

It is easy to see that the above set is convex and compact. Below we will 
study its extremal points. Let us shortly recall some necessary notations. Let 
A C X be convex set and X be some vector space. A point x € A is called 
extremal, if from 2x = x\ + X2, where x±, X2 € A and x\ ^ X2, it follows 
that x = x\ = X2- The set of all extremal points of a given set A is denoted 
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by extrA The Krein-Milman theorem asserts that any convex compact set 
on some topological vector space is a completion of the convex hull of its 
extremal points (see, for review [?]). 

Theorem 3.1. If A = (a^) <E extr U\ then an = V 1, a^ = V \ V 1 
Here and henceforth c = a V b means that c is either a or b. 

Proof. We prove this by induction with respect to the order of the set Ui. 
Let m = 2. 

First we prove that an = V 1 and 022 = V 1. Let < an < 1, then 
from an + 2ai2 + a 2 2 = 2 we get < a\ 2 < 1. Let us consider the following 
matrices 

/ a u + 2e a 12 - s ^ 
\ au-e a 2 2 
A = f an - 2e a 12 + e 

2 \ a 12 + e a 22 
Since < an < 1 and < ai2 < 1, then one can choose e such that 
Ai, A 2 G Ui. Now we obtain that 2 A = A 1 + A 2 . Therefore A £ extrUi. 
The last contradiction shows that an = V 1. By this analogy one can prove 
that 022 = V 1 Since an = V 1 and a 22 = V 1 it follows from A 6 Ui 
that aij = V \ V 1. One can easily see that the number of such matrices is 
4. That is the following matrices: 

;?)(?o)(H)(M 

And all of them are extremal. Thus, for A G extrXJi it is necessary and 
sufficient to be an = V 1, a^- = V \ V 1. However for m > 3 it is not true 
at all. 

Let us suppose that the conditions of the theorem is valid for all matrices 
of order less than m and prove it for matrix of order m. 

First we prove that if A € extr\J\ then an = V 1. Let us assume that 
it is not so, that is there exists some diagonal entry, which neither nor 1. 
Without loss of generality we may assume that < an < 1. 

We call a set a G / is said to be saturated, if Yl a ij = M- Corresponding 

i,jea 

minor to the saturated set we call saturated minor. For instance, I itself is 
saturated, but it is a proper set. Further, we mean only non proper saturated 
sets. 

If there is a saturated set, such that 1 G a then by the assumption of the 
induction we infer that an = V 1. Now remains the case when 1 ^ a for 
any saturated set a. 

Now we will prove that if < an < 1 then ajj = V 1 for all j. Indeed, 
if there is jo such that < a,j j < 1, then we put 

A' = {a'n = an +e, a'j j = dj jo ~ £ -> a ij = a %j f or a ^ other values of i,j} 
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and 



A" = {a'd = an — e, a" j = aj j +e, a'-j = aij for all other values of i,j} 

If jo is contained in some saturated set, then from the assumption of the 
induction it directly follows that a,j j = 0V1. If it is not so, then A', A" € Ui 
and consequently 2A = A' + A", that is A £ Ui, which is a contradiction. 

So if < an < L then ajj = V 1 for all j. Furthermore, since / is a 
saturated set, then there is ai j 7^ V \ V l(io 7^ jo)- Now we put 

A' = {a'n = on + 2e, j = aj j — e, aj,- = a^- for all other values of i,j} 



A" = {a'd = an — 2e, ^' j = a jojo+ e > a i'i = a ijf or a ^ other values of i,j} 

If {io,jo} is not contained in any saturated set, then A', A" G Ui and 
2A = A' + A". 

If {iojJo} is contained in some saturated set. Then by the assumption 
of the induction we get ai j = V \ V \{io 7^ jo), which contradicts to 

Oioj ^OVivl(io^jo)- 

Thus, all cases come to a contradiction. Therefore, an = V 1. Now we 
have to prove that if A € extrUi then a^ = V \ V \(i 7^ j). 

Let ai j 7^ V \ V \{%o 7^ jo)- If there is a saturated set a such that 
{*0)io} C a then by the assumption it follows that a^- = 0V^Vl(i 7^ j). Since 
/ is saturated, one can choose {ii, Ji} such that a{ 1 j 1 7^ V \ V l(zo 7^ Jo)- 
Also it follows that {h,ji} is not contained in any saturated set. 

Let us put 

A> = i a 'iojo = a io3o+ £ i a 'hji = a hji- £ , a ij = a ij f° r al1 other values of i,j} 



A " = i a i o j = a ioj - £ , a hn = ctiiji+e, a"j = aijfor all other values of i,j} 
The sets {io,jo} and {h,ji} is not contained in any saturated set. There- 



Corollary 3.2. Ifm = 3, then A <G extr U\ if and only if an = V 1, a^ = 



Moreover, \extrU\\ = 25 

Corollary 3.3. If A £ extr U\ then either A is stochastic or has a saturated 
minor of order m — 1. 



and 



and 



fore A', A" £ Ui and 2 A = A' + A". 



□ 



V \ V 1 and A / M. #ere M 
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Proof. Let A € extrXJi. Let us assume that the matrix A has no saturated 
minors of order m — 1. Then we have Y a ij < M f° r an y a such that 

m 

\a\ = m — 1. By taking into account Oy = m, we obtain 

Oji + ^2 > 1. 

By theorem it follows that a& = V 1, a^- = V \ V 1. Therefore an 
expression an + ^ a^ is an integer. Therefore 

an + ^ a^j > 2. 

By taking summation in the above inequality from i = 1 to i = m we find 

m mm 

y^(«M + ^ <%■) = 2 - > 2m - 

i=l j^i «J=1 i=l 

Hence 

m 

^2 aii - °- 

i=i 

and a^i = for all i = l,m. Further, an + Y a %j > 2 an d «ii = imply 



~Y aij > 1 for all i = 1, m. 

j'=i 

Since Y a ij = m then Y a ij = 1- The last means that A is stochastic. 
i,j=l j=l 

□ 

Theorem 3.4. Xei A = (aij)— be a symmetric and nonnegative matrix. 
For the existence of a stochastic matrix T = (tij) satisfying the equation (6) 
necessary and sufficient condition is 

ai i - H 

i,jea 

for all a C / = {1, 2, • • • , m}. 

Proof. Necessity. If T = (tij) is a stochastic matrix and (6) is satisfied, 
then 

Sufficiency. Let A G Ui . We have to show that existence of a stochastic 
matrix for which (6) is satisfied. 

First, we prove it for extremal points of Ui. Now let A £ extr\J\. Then 
from corollary 3.3 it follows that either A is stochastic or has a saturated 
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minor of order m — 1. If A is stochastic, then we take T = A since A is 
symmetric then (6) is satisfied. 

Now, let us suppose that A is not stochastic, then it follows that it has 
a saturated minor order m — 1. This case we prove by induction relatively 
order of the matrix. Let m = 2. In this case we have the following matrices 

?)■*-(! s)-*-(H)-*-(n> 

For Ai or Ai we take T = A\ and T = A2, respectively. 
For ^3 or A4 we take 

T= (i o) orT= (o 1)' res P ectivel y- 

Let us assume that assumption of the theorem is true for all matrices of 
order less that m. 

Since A has a saturated set of order m — 1, without loss of generality we 
may assume that this minor is j = j m _± . From the assumption of the 
induction there is a stochastic matrix for this minor which satisfies (6). Let 
it be V = (vij)ij = l,m-l. 

m—1 m m—1 

Since Y a-ij = m — 1 and Y a ij = mit follows that a mm + 2 Y a im = 

i,j=l i,j=l »=1 

1. From theorem it follows either a mm = 1 and aj m = for i = 1, m — 1 or 
aii = 0, ai om = \ ai m = for i = 1, m - 1 \ {io} 
For the first case we define T by 

-^{(^ij) • tij = Viji i) J ' = 1) Tn, 1 t m m = lj ^im = ^mi = 0> 1 = 1,771 — 1} 

and it is easy to see that above matrix is stochastic and satisfies (6). 
For the second case we define T by 

' ^ij — ^iji — 1,771 1 t m «o — 1' — ^™ — 0, 1 — 1,777,} 

and it is easy to see that above matrix is stochastic and satisfies (6). So, 
for extremal points of Ui the theorem has been proved. 

Now let A £ Ui From Krein-Milman(see [?]) it follows that A can be 
produced as the convex hull of its extremal points. 

Let 

k 

A = J2\ k A k . 

i=l 

k 

Here < A < 1, ^A^ = l and Ai G extrJJi, % = 1, k. 

i=i 

Let for matrices Ai correspond stochastic matrices Tj that satisfy (6). 
Then 

k 

T = ^2 ^kT k 
1=1 

is stochastic and satisfies (6). 
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□ 

Remark. It should be mentioned that stochastic matrix T (which is a 
solution of (6)) exists but not unique. 
Example. For the matrix 

/ 0,1 0,3 0,4 
,4= 0,3 0,1 0,5 
\ 0,4 0,5 0,4 

the set of the solutions of A = |(T+T') in the class of stochastic matrices 

is 

/ 0,1 a 0,9 -a 

T a = 0, 6 - a 0, 1 0, 3 + a 
\ a -0,1 0,7 -a 0,4 

here a G [0,1; 0,6]. 

Recall that a set A in R k is said to be a polytope, if it is nonempty, 
bounded and consists of intersection of a finite number of semispaces. 

Proposition 3.5. The set of all solutions of the equation 

A = \{T + T>) 
with respect to T forms a polytope. 

Proof. Let us rewrite the solutions of the equation in the following form: 

< % < 1, 

t%j ~\~ tji ^ t{j -\- tji ^ dij 

m m 

3=1 3=1 

From above given relations one can see that the set of all solutions of the 
equation is nonempty, bounded and it is intersection of semispaces. 

□ 

From (6) it follows that (Ax,x) = (Tx,x), here (•, •) is the inner product 
in R m . Indeed 

(Ax, x) = - ((T + T')x, x)=- [{Tx, x) + (T'x, x)] = (Tx, x). (7) 

Now let's return to the study of d.s.q.o. According to the theorem 2.1., 
if V : S™- 1 -» S™- 1 - d.s.q.o., then 

^ Pij,k < \ a \, Va G J, k = l,m 

moreover, for a = I an equality is held. 
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By putting 

Ak = {Pij,k)i,j=\~mi 

and using the theorem one can find a stochastic matrix such that 



Ak = \{T k + T' k ), k=l,m. 



From (7) one gets 



y ] Pij,k x i x j — (AfcX, x) — (TkX,x). 

Hence, if V is a d.s.q.o., then there are stochastic matrices T\,T 2 , ■ ■ ■ ,T n 
such that 

Vx = ((Tix,x), (T 2 x,x), ■■■ , (T m x,x)) . 
Now we prove the following inequality: 



m— 1 



min Xi < (Tx,x) < max Xj, Vx G S 

l<i<m l<i<m 

or in accepted notations: 

x [m] < (Tx,x)< X[1] , VxeS m -\ (8) 
here T— is a stochastic matrix. 

m 

Indeed, if tij > and Y Uj = 1) then 

m 

X[m] UjXj < X[l] 

3=1 

for all x € R m . In particular,for x € S*™ -1 we get 



m 

(Tx,X^) — ^ ^ tijXiXj — ^ ^ x i ^ ^ ^ij^j' 
i,j=l i=l j=l 

m 

By Xj > 0, Y x i = 1 an d above inequality we get the required inequality 

i=l 

(8). 

Theorem 3.6. Let A = (aij) be a nonnegative symmetric matrix. Then for 
the fulfillment of the inequality 

X[ m ] < (Ax,x) < X[i] Vx e 5 m_1 , 

it is necessary and sufficient that 

aij < \a\ Va € /, fe = l,m 

moreover, for a = I an equality is held. 
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Proof. Sufficiency straightforwardly follows from Theorem 3.4 and the in- 
equality (8). 

Necessity. Consider a set a C / and put x° = (x\, x®, ■ ■ ■ , £ m ), with 

x\ = - — -, i G a, = 0, i a. 
\a\ 

Then for \a\ < m we get X[ m ] = 0, x^j = j^. Therefore, one has 

m 1 1 

< (Ax, X) = OijXiXj = ^2 dij < 

hence 



\ct\ a 



If \a\ = m, then X[ m ] = ^[i] = ^- Consequently 

11^ 1 

— < — o V a»i < —j 



m L — ' m 



rn 



that is Yl a ij = m - O 
By Theorems 3.4 and 3.5 one can deduce that an inequality 

x [m] < (Ax,x) < x {1] Vx G 5 m_1 (9) 

here A = (aij) is a nonnegative matrix, is equivalent to the existence of 
stochastic matrix T = (Uj), such that 

A + A' = T + T'. 

Since we are going to study sufficient conditions for doubly stochasticity 
of q.s.o. from the proved theorem one arises the following 

Problem 3.7. For which symmetric matrices A the following inequality 
holds 

X[ m ] + x [m-l] H 1" x [m-k+l] < { Ax , x ) < x [l] H ^ x [k] , (10) 

/or a// x G S' m ~ 1 , /tere k— natural number 1 < k < ml 

For = 1 an answer is given by theorem 3.3. To solve the problem we 
give some notations: 



T k = {T = (Uj), i,j = l,m:Q< % < 1,^2 Uj = k}, 1 < k < m. 

3=1 

U fc = {A = (aij) : aij = a jU A=^(T + T'), T G T k }, 1 < k < m. 
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The set Ufc is called a symmetrization of TV Evidently, Ti— the set of 
all stochastic matrices. 

Theorem 3.8. For A € the inequality (10) holds. 

Proof. Let x^ = (xm,--- , xr m i) — be a nonincreasing rearrangement of x 

m 

and Ai, • ■ ■ , X m be an arbitrary numbers such that < Aj < 1, ^ Aj = k. 
Consider the following sum 

AiX[!] + A 2 X[ 2 ] H h A m X[ m ]. (11) 

Let, i < j. Replacement of the coefficients Ai, ■ • ■ , Aj, • • • , Aj, • • • , \ m into 
Ai, Aj + £,•■■, Aj — £,•■■ , A m we will call backward shift, if e > and 
forward shift if e < 0. Shift is called admissible, if it preserves the condition 
< A£, here A^ are the coefficients, obtained in the result of shift. 

Since x^ > xyj, then it is clear that under the admissible backward shifts, 
the sum (11) does not decrease, and respectively, does not increase under 
the admissible forward shifts. 

It is easy to see that an admissible backward shift is possible up to ob- 
taining the collection 1, 1, • • • , 1, 0, 0, • • • ,0, and to the right up to obtaining 
the collection 0, 0, • • • , 0, 1, 1, • • • ,1, here in both cases the number of ones 
is k. That's why 

x [m-k+l]+ x [m-k+2]-\ \-X[m\ < AiXjx] +A 2 X[ 2 ] H hA m X[ m ] < X[!]H |-X[ fe ] 

(12) 

m 

for an arbitrary x <G R m and < Aj < 1, ^ Aj = A;. 

i=i 

Now let x e S m ~ l and A <G Ufc. 
Choose T € Tfc,such that 

A = \{T + T') 

Then 

m mm 

(Tx, x) — ^ ^ tijXiXj — Xj(^^ ijjXj) 
j,j=i i=i j=i 

m 

Since < ijj < 1 and ^ tjj = k, from (12) one gets the following 

m m fc 

i=m— fc+1 j=l i=l 

m. 

According to Xj > 0, ^ Xj = 1 we obtain 
i=i 

m A; 
i=m— fe+1 i=l 
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or 

m k 

x [{l <(Ax,x)<^2x li] 

i=m—k+l i=l 

So if A G U fe then (10) holds. □ 

Let B be the set of all d.s.q.o. By putting Ak = {Pij,k}(i j=Tm) we rewr it e 
an operator V in the following form 

V = {A l \A 2 \---\A m ). (13) 
Then theorems 2.1, 3.4, 3.5 and 3.6 imply that for VeB the conditions 

A k G Ui, k = l,m (14) 
are necessary and the conditions 

Vac/, J> fe eU H (15) 

fcSa 

are sufficient. 

Whether the conditions (14) is sufficient to be V G B? We give a positive 
answer for the question in next section. 

Now we show some simple properties of the which will be helpful. Let 
E be the matrix of m x m with all entries equal to unit. 

Theorem 3.9. The following assertions hold 

i)A(£U k ^E-A£ U m ^ k ; 
ii) U k n[/i = 8, k + I; 
Hi) U rn = {E}; 
iv) Ae U k ^j-Ae Up, 1 < p < k; 
v) U k + UiD U k+i , k + l<m. 
Proof, i) Let 1 < k < m — 1 and A G U k . From 

A= l -(T + T>) 

m 

here T = (Uj), < < 1, Yl Uj = k one gets 

E - A = ^[(E - T) + (E - T)'] 

It is obvious that E — T G T m - k , therefore E — A G U m _fe. 

ii) It follows from the definition of TJ k . 

iii) Since A G U& then < a^- < 1 for all k = l,m. If A G U m then 

TO 

Y aij = m 2 and therefore aij = 1 that is U m = 



14 



RASUL GANIKHODZHAEV AND FARRUH SHAHIDI 



iv) Let A G U k and A = \{T + T'). Then 

moreover |T G T p . Consequently, |A G U p . 

v) Let Hl<m,ie U fe+ i and A = \{T + T') By denoting 

A 1 = - k - (T + T'), A 2 = — 1 J T + T') 

we obtain A = A\ + A2, furthermore A\ G Ufc, A 2 £ Uj. The last implies 

U t + UjDU w , k + l<m. 

It should be mentioned that the last inclusion is strict for all k, I > 1 Indeed, 
let / be the identity matrix. Clearly that / G Ui, but I + I £ XJ 2 - □ 

4. Extreme point of the set of doubly stochastic operators 

and blrkhoff's problem 

According to the classic result of Birkhoff (see [2]) extreme points of the 
set of doubly stochastic matrices are permutation matrices, i.e. matrices 
having exactly one unit entry in each row and exactly one unit entry in 
each column, all other entries being equal to zero. It is interesting to know 
an answer for the similar problem about the set of all doubly stochastic 
nonlinear operators. Therefore, to investigate such a problem it is better 
first to start with the set of all d.s.q.o.'s, since such a set contains as a 
subset of the set of all doubly stochastic matrices. In this section we are 
going to describe extreme points of the set of d.s.q.o. At first, we prove 
necessary and sufficient conditions for q.s.o. to be d.s.q.o. (it was proved in 
the last section but in this section we show that necessary and sufficient will 
coincide). 

The set of all d.s.q.o. we denoted by B. Directly from definition it follows 
that 



V = (Ai\A 2 \ ■ ■ ■ \A m ) G B 4» V = (^(1)1^(2)1 • • • K (m)) G B 
here tt is an arbitrary permutation of the index set / = {1, 2, • • • , m} 

Theorem 4.1. Assume that A± G U±. Then one can choose A 2 , • • • , A m G 
t/i such that V = (A 1 \A 2 \ ■ ■ ■ \ A m ) G B. 



Proof. Take 



A 2 = A 3 = ■ ■ ■ = A m = —}—(E - Ai) 

m — 1 



here E is the matrix with all entries equal to unit. According to the theorem 
3.8 we have 

E-A 1 € U m _i 
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and 

—L-(E-Ai) eUi 

m — 1 

Let a stochastic matrix Ti G Ti be a solution of the equation 

A x = \{T+r) 

Now we put 

T 2 = T 3 = • • • = T rn = — !—(£ - Ti) 

m — 1 

we show that the sum of any k of the matrices Ti, T2, • • • T m belongs to 
TV Indeed, let T\ = {ty} such that < < 1. Then 

T2 = T 3 = • • • = T m = -{1 — tij} 

m — 1 

From < Uj < 1 one gets 

0<— ^-(1-%)<1 (16) 

and 

k — 1 

0<% + -(1 - tij) < 1. (17) 

If the sum contains T\ then our assertion follows from the inequality 
(17), otherwise from the inequality (16). So the sum of any k matrices of 
Ti, ■ ■ ■ ,T m belongs to T k . 

Now we return to Ai, ■ ■ ■ ,A rn and obtain that the sum of any of them 
belongs to Ufc. Therefore we get that 

V = (A 1 \A 2 \---\A m )eB 

□ 

v 

It is clear that if Ti,-- - ,T p G Ti(p < m) and J2 T i G t p then for 

i=l 

any k < p the sum of any k matrices of Ti,-- - ,T p belongs to T^. This 

p 

implies that from A±, ■ ■ ■ , A p G Ui, ^2 Ai <=L JJ p it follows that the sum of 

i=i 

any k(k < p) matrices belongs to Ufc. Using the last assertion we get the 
following 

p 

Corollary 4.2. If A\,--- ,A p G Ui, ^ Ai ^ U p (p < m) then one can 

i=i 

choose A p+ \, ■ ■ ■ A m G U\ such that 

V = (A 1 \A 2 \---\A m )eB 
So conditions (15) can be changed to conditions: 

m 

ijGUi, i = T~m, J2 A i = E ( 18 ) 
i=i 

since U m is E. Therefore 
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Corollary 4.3. The conditions (18) are necessary and sufficient for q.s.o. 
to be d.s.q.o. 

Corollary 4.4. The set of all d.s.q.o. forms a convex poly tope. 

Proof. We recall that finite intersection of nonempty, bounded and closed 
semispaces is called convex polytope (see [6]). For any q.s.o. there is a cubic 
matrix which can be embedded in the space R m . One can embed it this 

rrp" (m — 1) 

cubic matrix in R 2 . Furthermore, each of the condition (18) defines 

mP" (m— 1) 

closed semispace in R 2 . If we denote the set of all d.s.q.o. by B, from 
the corollary 4.3 it follows that consists of intersection of these semispaces. 

vr? (rn — 1) 

It is clear that this intersection is nonempty and bounded set in R 2 , 

m 

since it lies in positive ortant and in hyperplane Pij,k = m2 ■ Due to 

i,j,k=i 

Grunbaum [6] we conclude that B is convex polytope. □ 

In [3] it was conjectured that an operator V = (A1IA2I • • • \A m ) is extremal 
point of B if and only if A; L G extrU, Vi = 1, m. We show that in general it 
is not true. 

Theorem 4.5. Let V = (Ai\--- \A m ) G extrB. ThenV w = (Ar(i)l ■ ■ ■ |Ar(m)) G 
extrU for any permutation tt of the index set {1, 2, • • • , m}. 

Proof. Let V n <£ extrBThen 3V',V" G B, V ^ V", such that 2V n = 

V + V" . Let V = (A[\ ■ ■ ■ \A'J, V" = (A'{\ ■ ■ ■ Then (2A 7r(1) - A[ - 
A'{\ ■ ■ ■ |2A^( m) - A' m - A'jj = 0. Since the matrices 2A W ^ - A \ - A'( are 
symmetric, then 2A n ^ = A\ + A", therefore 2Ai = A' 7r _ 1 ^ + A'^_ l( ^y Since 

m m 

i=i i=i 

then 

m m 

Denote W = {A'^J- ■ ■ |<_ 1(m) ) and W" = ■ • • |<_ 1(m) ). Then 

we get W, W" G B, W ^ W" 2V = W + W" which contradicts to 

V G" extrB. 

□ 

Theorem 4.6. Let V = (Ai\ • • • \A m ). If for any permutation ir of any m—1 
elements of {1, 2, • • • m} we have A w ^ G extr U , then V G extrB. 

Proof. Using the above theorem it is enough to show that if A±, A2, ■ ■ ■ , A m -\ G 
extrXJ, then V G extrB. Let us assume that V £ extrB. ThenBV' ,V" G 
B, V / V" such that 2V = V + V'.Let V = (A[\--- \A' m ), V" = 
KI---KJ. Then 

(2 Al - A[ - A'{\ ■ ■ ■ \2A m - A' m - A" m ) = 0. 
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So 2Ai - A\- A'l = Vz = l,m. Since A; £ extrf/, Vz > 2, then A\ = 

m m 

A'/, Vi > 2. From the equality E ^ = E ^" = ^ we S et = A i- That's 

i=i i=i 

why V' = V", which contradicts to V £ extrB. □ 



Corollary 4.7. If Ai £ extrU, Mi = l,m, then V £ extrB. However, the 
converse case is not true 

Example : Consider an operator 

x' x = XlX 2 + X1X3 + x 2 x 3 
X2 = X2 + X3 + X1X3 
x 3 = x\ + XlX 2 + x 2 x 3 

Corresponding matrices have the following view 






It is easy to see that two of them are extremal. Then from the above 
theorem it follows that V £ extrB. 

Whether the conditions of the theorem are necessary and sufficient? In 
two dimensional simplex the problem is solved positively. 

Theorem 4.8. Let V = (A 1 \A 2 \A 3 ) £ B. V £ extrB if the only if at least 
2 of the matrices Ai,A2,A% are extremal of U±. 

Proof. For the proof of the given theorem we refer to the theorem 5 from 
[4], which says that if V £ extrB, then entries of the matrices A\,A2,Az 
either or \ or 1. Therefore, if V £ extrB, then either A\,A 2 , A3 £ extrU 
or A 1 ,A 2 ,A 3 = M. 

The case A\ = M, A 2 = M,A% = M is impossible because of A\ + A 2 + 
A 3 = E. 

It can be easily shown that if two of the matrices A±,A 2 , A3 is M then V 
is not extremal. Therefore we can conclude that at least two of the matrices 
A\ , A 2 , A3 are extremal. □ 

Corollary 4.9. For m = 3 we have |extrB| = 222. 

Proof. Let V £ extrB, V = (A1IA2IA3). From the theorem 4.8 it follows 
that A 2 ,Az £ extrU up to permutation. From A\ + A 2 + A3 = E one gets 
that either A\ £ extrJJ or A\ = M. 

Let A\ £ extrU. From corollary 3.2 we know all extreme points of U 
Therefore we can choose those triples of extreme points, sum of which is E. 
The number of such triples is 31. For the case of A\ = M, the number of 
such triples is 6. 

Consequently, |extrB| = 37 x 3! = 222. □ 
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